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We present an escape rate theory for current-induced chemical reactions. We use Keldysh nonequi- 
librium Green's functions to derive a Langevin equation for the reaction coordinate. Due to the 
out of equilibrium electronic degrees of freedom, the friction, noise, and effective temperature in 
the Langevin equation depend locally on the reaction coordinate. As an example, we consider the 
dissociation of diatomic molecules induced by the electronic current from a scanning tunnelling mi- 
croscope tip. In the resonant tunnelling regime, the molecular dissociation involves two processes 
which are intricately interconnected: a modification of the potential energy barrier and heating 
of the molecule. The decrease of the molecular barrier (i.e. the current induced catalytic reduc- 
tion of the barrier) accompanied by the appearance of the effective, reaction-coordinate-dependent 
temperature is an alternative mechanism for current-induced chemical reactions, which is distinctly 
different from the usual paradigm of pumping vibrational degrees of freedom. 



I. INTRODUCTION 

The recent advances in nano fabrication make it pos- 
sible to use the scanning tunnelling microscope (STM) 
as a "nanoscale chemical reactor". The tunnelling cur- 
rent from the STM tip can selectively break or form 
chemical bonds^"'^ and initiate chemical reactions of the 
reactantsi^i^ 

The rates of chemical reactions can be routinely com- 
puted for molecular systems in thermodynamic equi- 
librium. Usually one relies on the Born-Oppenheimer 
approximation^^^ and may include various nonadiabatic 
effects, when the electronic energy levels are not well sep- 
arated or coupled to the continuum of the environment 
statesi^rJ^ Let us now assume that the electronic system 
is driven out of equilibrium. As an example for an out 
of equilibrium system we consider a molecular junction 
(molecule attached to two macroscopic metal electrodes 
held at different chemical potentials) or a molecule ab- 
sorbed on a metal surface under a scanning tunnelling 
microscope tip. As electric current is flowing through the 
molecule, considerable amounts of energy are dissipated 
and partly passed from electronic to nuclear degrees of 
freedom — in linear response the total dissipated power 
is proportional to IV. Thus, since the current is in the 
range of nanoamperes and the voltage is a few volts, a 
significant energy (~ 0.2 Hartree per nanosecond) is dis- 
sipated in total and the energy which is dissipated in the 
molecule per nanosecond can be comparable to typical 
barriers for chemical reactions. 

Generally speaking, the electrons produce not only the 
standard adiabatic force but also give random fluctua- 
tions and viscosity to the nuclear dynamics..^^ Close to 
equilibrium, the latter are related by the fluctuation- 
dissipation theorem,A2. but far from equilibrium the noise 
is no longer balanced by the viscosity. The absence of the 
fluctuation-dissipation relation as a direct consequence of 



nonequilibrium opens new possibilities in chemistry, such 
as unusual ways to catalyse chemical reactions. 

In this paper, we develop a theory for current-induced 
chemical reactions. The main physical assumptions of 
our theory are that the electronic dynamics is much faster 
than the nuclear motion and that the reaction coordi- 
nate is a classical variable. This results in Markovian dy- 
namics of the reaction coordinate described by a Fokker- 
Planck equation. This Fokker-Planck equation has in- 
teresting new features. Namely, the nonequilibrium po- 
tential energy surface depends on the electronic current 
flow through the molecule and the effective temperature 
produced by the nonequilibrium electrons on nuclear de- 
grees of freedom is no longer constant but becomes a 
function of the reaction coordinate. In other words, the 
out of equilibrium electrons not only change the profile of 
the potential energy surface for the reaction coordinate 
but locally heat it with different rates. For the paradig- 
matic example of the dissociation of diatomic molecules, 
we show that current-induced chemical reactions involve 
two interconnected processes: a modification of the po- 
tential energy barrier and heating of the molecule. This 
new mechanism for current-induced chemical reactions, 
complements the familiar paradigm of pumping vibra- 
tional degrees of freedomJ^ii^ The principal difference 
between this paper and previous work^'^ lies in lifting the 
strong assumption that the temperature of the nuclear 
degrees of freedom is exactly the same as that of the 
equilibrium electrons in the metal surface. 

The remainder of the paper is organized as follows. In 
Sec. 2, we describe the Langevin equation for the reac- 
tion coordinate. Section 3 presents the calculation of the 
reaction rates via the solution of a Fokker-Planck equa- 
tion. Conclusions are given in Sec. 4. Some technical 
aspects are relegated to appendices. We use atomic units 
throughout the paper. 



II. LANGEVIN EQUATION FOR THE 
REACTION COORDINATE 

As a model system, we consider a diatomic molecule 
attached to two metal electrodes, say, a metal surface on 
one side and an STM tip on the other. The molecule is 
modeled by one electronic spin-degenerate molecular or- 
bital with single particle energy e{x) + Vg, which depends 
on the bond length x and the gate voltage Vg . The bond 
length X can be considered as the reaction coordinate. 
The Coulomb interaction between electrons is accounted 
for by a charging energy Uc (x) which is a function of the 
bond length x. The reaction coordinate x is considered to 
be a classical variable with corresponding momentum p 
and a reduced mass m (taken as that of the H2 molecule, 
m = 918 a.u.). The nuclear Coulomb repulsion energy is 
Vn{x) — 1/x. Then the molecular Hamiltonian is 

1 

Hm = {e{x)+Vg)22"-t^'^+2^c{x)n-^ni + — + VNix). 

a 

(1) 

Here a\{aa) creates (annihilates) an electron with 
spin a =t, i in the molecular orbital and Ua = al-Qa- The 
details of the parametrization of the molecular Hamilto- 
nian are given in appendix A. 

The complete Hamiltonian of the molecular junction 
consists of the molecular Hamiltonian ([IJ , the Hamilto- 
nians for noninteracting left and right electrodes, and the 
molecule-electrode interaction: 

H = Hm + ^ ekal^.acrk + ^ {tkal^aa + h.c), 

a.keL,R. a,keL.R 

(2) 

where a^j,(acrfc) creates (annihilates) an electron in the 
single-particle state ak of either the left (L) or the right 
(i?) electrode. The electron creation and annihilation op- 
erators satisfy standard fermionic anticommutation rela- 
tions. The electron occupation numbers in both elec- 
trodes obey the Fermi-Dirac distribution fL,R{£k) — 
[1 _l_ g(efc-/'i,.H)/Tj-i^ The chemical potentials in the elec- 
trodes, and /Zfl, are assumed to be biased by the ex- 
ternal symmetrically applied voltage V = fiL — fJ-R, which 
we take to be positive, and hl,r = ±0.5V. The Fermi 
energies of the electrodes are set to zero. We also assume 
that the tunnelling matrix element tk is spin indepen- 
dent. 

The coupling with the electrodes broadens the molec- 
ular level, with the width given by the imaginary part of 
the electrode self-energy, 

r(w) = Fl + F;?, = -Im J2 \tk\ViuJ-ek+tO+). (3) 

keL,R 

In what follows, for the sake of simplicity, we will use 
the wide-band approximation for the electrodes, i.e., the 
imaginary part F of the self-energy is energy independent 
and the real part vanishes. The total width of the molec- 
ular level is fixed in our numerical calculations, F = 0.05, 
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FIG. 1: Upper part: Possible experimental setup. The STM 
tip is positioned above a hydrogen molecule (represented as a 
sphere). An electron coming from the tip onto the molecule 
then tunnel into the conducting surface. A gate electrode is 
used to tune the molecular energy level(s) and to control the 
current flow between the STM tip and the conducting surface. 
Lower part: Schematic illustration of the model. By varying 
the gate voltage the molecular energy level can be shifted 
near the chemical potential of the electrodes. The molecular 
orbital energy £(2;) and electron-electron interaction energy 
Uc(x) depend on the bond length x. 



but we vary the relative contributions of left and right 
electrodes. 

We neglect the effect of the electric field on the 
molecule - for diatomic molecules lying fiat on the metal 
surface the electric field is perpendicular to the bond 
length and does not infiuence the dynamics. In other 
cases, the electric field may give an additive contribu- 
tion which can be easily included into the model for par- 
ticular molecule-electrode geometries and inter-electrode 
distances. 

We begin with the Ehrenfest coupled electron-nuclei 
dynamical equations 

ihp{x,t)^[H,p{x,t)] (4) 
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V = 0.01 ^-"'"''^ 


V = 0.05 -*»-»***^ 




FIG. 2: Effective single particle energy level Ecff (a^) and pop- 
ulation n{x) as a function of the reaction coordinate x for 
different values of the asymmetry coefficient g — Tl/Tr 
and applied voltage V. We chose eefr(2:min) = —0.36, where 
Xmin = 1-66 corresponds to the minimum of the equilibrium 
potential energy of the molecule. The Fermi energies of the 
electrodes are set to zero. 



the Hamiltonian ([T]), i.e., we replace it by 
H'^/ = £eff(x) ^ + Vn{x) + 



where 



e{x) + Vg + -Ucix)n{x) 



(8) 



(9) 



is the effective single particle energy and n{x) — 
TT[p{x)na-] is the nonequilibrium population of the elec- 
tronic level of the molecule. This can be computed by 
means of the NEGF, see Eq. (jB17|) in Appendix B. In 
Appendix B we also present the details of the derivation 
of the explicit expressions for ^{x) and D{x). 

For low temperatures, T — > 0, the Fermi-Dirac electron 
distributions in the electrodes fL,Riy>) can be approxi- 
mated by a step like function, and the level population, 
friction coefficient, and noise intensity can be written as 



+ Tr arctan ('fsii^Lj^^) i> (iq) 



2r 

= — {TLQ{^^L,x) + TrQ{pr, x)} , (11) 



dx 



Tr[p{x,t)H] 



(5) 



We eliminate the electronic density matrix from this 
equation assuming that the electronic degrees of freedom 
vary much faster than the nuclear motion. Physically 
this means that the oscillation period T of the reaction 
coordinate is much larger than the tunneling time for the 
electrons, i.e., T ^ l/F. Below we will see that this as- 
sumption is fulfilled for the model under consideration. 
This eliminates the direct time dependence from the elec- 
tronic density matrix p{x,t) ^ p{x), which now depends 
parametrically on time through the reaction coordinate. 
The result is the following Langevin equation for the re- 
action coordinate, 



-Tr[p{2 



dH 

' dx^ 



^{x)x + Sf{x,t). 



(6) 



Here p{x) is the nonequilibrium density matrix, ^(x) is 
the frictional force (viscosity), and Sf{x,t) is the random 
force (noise) taken in the Gaussian form 

{Sf{x, t)) = 0, {Sf{x, t)Sf{x, t')) = D{x)5{t - 1'), (7) 

The derivation of the Langevin equation via Keldysh 
nonequilibrium Green's functions (NEGF) is presented 
in Appendix B. 

To obtain expressions for the frictional force and noise 
intensity D{x) we apply the mean- field approximation to 



D{x) — / dujQ{uj,x) 



where the function Qiuj^x) is defined as 



Q{uj,x) 



(c^-eeff(x))2+F2 



(12) 



(13) 



Note that in the derivation of Eq. pT|) we have utilized 
{-duifa) = fa{l - fa)lT and lim {-d^fa) = S{oj - p.^). 

Equations (0 HH [HI ITU)) are the main equations of our 
model. Their solution provides us with the parameters 
of the Langevin equation ([6]). Note that ees{x) in (flO)) 
itself depends on n{x), through Eq. Therefore, both 
equations should be solved self-consistently. The main 
variables in our model are the applied voltage V and the 
asymmetry coefficient g = Tl/Tr. The asymmetry coef- 
ficient g controls the relative strength of the coupling to 
left and right electrodes. Replacing ghy l/g is equivalent 
to reversing the applied voltage. 

In Fig. [2] the effective single particle energy ecs{x) and 
the population n(x) are depicted for two values of the ap- 
plied voltage, small {V = 0.01) and large {V — 0.05), and 
for different values of the asymmetry coefficient g. The 
reference point for the effective molecular orbital energy 
is chosen to be eeff(a;i„in) — —0.36 (relative to the elec- 
trode Fermi energy), where Xmin = 1-66 determines the 
minimum of the H2 ground state energy [see Ea. (jA3p ]. 
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FIG. 3: Friction coefficient E,(x) and noise intensity D(x) as 
a function of the reaction coordinate x. All the parameters 
are the same as those in Fig. (2] 



This reference point can be shifted by the real part of 
the electrode self-energy, if we revoke the wide-band ap- 
proximation, or by the application of a gate voltage Vg. 
The dependence of the results on gate voltage will be 
discussed below. As seen from the figure, at small values 
of the reaction coordinate, i.e., when the electron level 
is well below the electrode Fermi levels, both £cff(a;) and 
n(x) are nearly independent of the applied voltage and 
the asymmetry coefficient. The situation changes with 
increasing x, i.e., when £off (a;) comes into resonance with 
the electrode Fermi levels. In that case, the electron pop- 
ulation is strongly affected by the asymmetry coefficient 
and the effect becomes more pronounced at larger values 
of V . As a result, at large x and V the effective single 
particle energy Eeff (a;) grows with increasing g. 

In Fig. [3] the friction coefficient ^{x) and the noise in- 
tensity D{x) are plotted for the same parameters as those 
in Fig. [5] The friction coefficient grows with increasing 
X and reaches its maximum value when Easix) is close 
to the chemical potentials of the electrodes, i.e., when 
the resonant regime is approached. The friction ^(a;) 
becomes smaller when the molecule starts to fall apart 
(due to the decrease of e^g(a;)). As seen from the figure, 
the friction coefficient grows as the asymmetry coefficient 
is reduced and this dependence is more pronounced for 
larger values of V . Besides, for small values of g the fric- 
tion is more sensitive to changes in the applied voltage. 

The noise intensity T}(x) also reaches its maximum 
value when ecfi(a;) is close to the chemical potentials 
of the electrodes. However, in contrast to the previ- 
ous case, D(x) is maximal for symmetric coupling to the 



FIG. 4: Upper panels: Effective temperature Tcff as a func- 
tion of the reaction coordinate x for two values of the applied 
voltage X. Lower panels: Effective temperature T^a as a func- 
tion of the applied voltage V for two values of the reaction 
coordinate x. All the parameters are the same as those in 
Fig.El 



electrodes, i.e., when g = 1. Moreover, because of the 
integral in Eq. (fT2|) the order of magnitude of D{x) is 
proportional to the applied voltage. 

Motivated by the fiuctuation-dissipation theorem, we 
can also introduce an effective "temperature" which de- 
pends on the reaction coordinate. 



(14) 



2'?(a;) YLQ{^^L■,x)+YRQ(^[lR,x)' 

const for fiL > io > fiR we can 



Assuming that Q{uj,x) 
approximate (1141) by 
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(15) 



The effective temperature ([T5|) has a maximum of V/A 
when 5 = 1. Note that the current is also maximal when 
the coupling to the electrodes is symmetric. In the upper 
panels of Fig.|4]the effective temperature Tcff (a:) is shown 
for the same parameters as those in Fig. [5] for two values 
of the applied voltage. As seen, for small applied voltages 
rcff(x) is nearly independent of x and can be approxi- 
mated by Eq. (fT5|) . When the applied voltage increases 
the effective temperature also grows and Tes becomes de- 
pendent on the reaction coordinate in a rather nontrivial 
way. The lower panels of Fig. |4] depict the effective tem- 
perature as a function of V for two particular values of 
the reaction coordinate. The point x = 2.0 is close to the 
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FIG. 5: EfTective nonequilibrium potential L'cff(2;) as a func- 
tion of the reaction coordinate x for two values of the applied 
voltage V. The panels correspond to different positions of 
the effective molecular orbital energy relative to the electrode 
self-energy £off(a;min) ~ —0.36 + Vg. All other parameters are 
the same as those in Fig. [2] 



minimum of the effective potential, while x — 5.0 is near 
the top of the barrier. We emphasize that by reversing 
the applied voltage polarity (i.e., replacing g hy 1/g) we 
can vary the effective temperature and the effect is more 
pronounced the larger the value of V. 



III. CALCULATIONS OF THE REACTION 
RATES 

A. Fokker-Planck equation for the reaction 
coordinate 



Since dH/dx 



'cffV 



c) X^CT '^o" + ^Ar(^) '^^ mean- field 



approximation, Eq. ^ becomes 
mi = -2e'^s{x)n{x) - V^{x) ~ ^{x)x + 6f{x, t). (16) 



The Langevin equation (jl6p is equivalent to the Fokker- 
Planck equation for the distribution function F{x,p,t), 



m dp 



2 9p2 



(17) 
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TABLE I: The period of motion (Eq. ([1^1) ) in the nonequi- 
librium effective potential {V = 0.05) for different values of 
the asymmetry coefficient g and the gate voltage Vg. The pa- 
rameter k defines the energy, k — {E — Emin) / (Smax — i^min)- 
All other parameters are the same as those in Fig. (2] 



Here the effective nonequilibrium potential energy sur- 
face is defined via integration of the nonequilibrium force 
in the Langevin equation. 



U,f[ix) = Vn{x) + 2 / e'^^{y)n{y)dy 



(18) 



In Fig. [5] the nonequilibrium effective potential is 
shown as a function of the reaction coordinate for differ- 
ent values of the asymmetry coefficient applied voltage 
bias V, and gate voltage Vg. The current flow through the 
molecule changes the height of the potential barrier: for 
g < 1 the current reduces the barrier height and for g > 1 
the current increases it. The explanation is the follow- 
ing: when the coupling to the left electrode is stronger 
{g > \) the current pumps electrons into the molecule 
and enhances the chemical bond. Conversely, when the 
coupling to the right electrode is stronger the current 
depletes the molecular electrons, thereby weakening the 
chemical bond. Thus, by reversing the voltage polarity 
we can vary the height of the potential barrier. But it 
should be stressed that the effective temperature T^s also 
changes with reversing the voltage polarity and the bar- 
rier growth is accompanied by a temperature increase. 
We also see from Fig. [5] that the height of the potential 
barrier can be decreased by increasing the gate voltage. 
The resonant tunneling regime Vg = and Vg — 0.05 
is physically most interesting, since the nonequilibrium 
potential energy profile exhibits a clear barrier between 
product and reactant states. We focus on the resonant 
tunnelling regime in our calculations of the reaction rates. 

Now we want to reduce the Fokker-Planck equa- 
tion p?)) to the Smoluchowski equation for the distri- 
bution function which depends either on the reaction 
coordinate x (overdamped motion) or on the energy 



~ 2m~^ t/cff (a;) (underdamped motion). To find which 
case is realized for the considered model we compare the 
oscillation period along the reaction coordinate and the 
relaxation time due to friction. The oscillation period is 
given by 



, , n dx 
T{E) = 2 1 — = 2 



X 



dx 



-1 ^^{E~U,^{x)) 



(19) 
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where Xi and X2 are left and right turning points, i.e., 
E = Ucsixi,2)- If ~ f7eff(2;min) then T « 2tt/loo, where 
cjo = \/U^'g{xniin)/m is the oscillation frequency near the 
bottom of the effective potential In Table U the oscilla- 
tion period is computed for various values of the bias 
voltage, gate voltage, and asymmetry coefficient. Since 
T ^ 1/r (r = 0.05) the assumption behind the deriva- 
tion of the Langevin equation (appendix B), namely the 
assumption that the electronic degrees of freedom are 
much faster than the motion along x is fulfilled. In Fig. [3] 
the friction ^(x) is shown for different values of asymme- 
try coefficient g and applied voltage V. It is evident 
that the relaxation time due to the friction m/^(x) is 
much larger than the oscillation period. Therefore the en- 
ergy dissipation per period of the motion is small, which 
means that we deal with underdamped motion. 



B. Calculations of the escape time 

The solution of the Fokker-Planck equation is not triv- 
ial in our case, since the effective temperature depends 
on the reaction coordinate. Using the method described 
by Coffey et alJ^ the Fokker-Planck equation is reduced 
to the equation for the distribution function for the en- 
ergy, P{E,t) 



where 



l_d_ 

2dE 



where 



and 



TiE) J,, 



D{E)-^P{E,T) 



^{x)xdx 



D{E) = 



TiE) J,, 



D{x)xdx. 



By introducing 



a{E)^\-^D{E)~m 



(20) 



(21) 



(22) 



(23) 



we re-write this equation as the Smoluchowski equation 
^P(E,t) = - A [„(i5)p(i5,i)]+l^ [D{E)P{E,T)] . 

The same equation was obtained in Ref. 1171 . though 
the theory was applied only to a nano-mechanical har- 
monic oscillator and thus did not consider dissociation 
processes. Then, the mean time for a particle with 
initial energy -Emin = Uaeixmin) to arrive at the top 
Emax — C/cff (2^max) of the potential barrier is given by 
(see, for example, Eq. XII. 3. 7 in Ref. IT8| ) 



2 / d£;e*(^) / 



dE' 

D{E')' 



HE') 



(25) 



^ ' Jeo D{E') 



(26) 



Note, that the choice of i?o is not relevant since it does 
not contribute to the escape time ([25]) . Using rela- 
tion (l23l) we can rewrite (1251) as 



dE 
D{E)' 



e-^^^'UE' 



where 



J En J-cS( E ) 



(27) 



(28) 



in terms of the effective temperature Ti:.s{E) = 
D{E)/2CiE). 

Equation ([571) gives us the exact escape time (mean 
first passage time) for the Brownian particle in the un- 
derdamped regime with coordinate-dependent effective 
temperature. In our particular case it can be further 
simplified, because our calculations show that the effec- 
tive temperature Tcs{E) depends only slightly on energy 



when Er, 



< E < E„ 



and can be approximated by 



Eq. (HH)., i.e., Tcs{E) = T{V,g). With this approxi- 
mation, the integral becomes ^'(E') w E/Tcs{E). 
Substituting this into Eq. ((27)) we obtain 



exp 



E — Emir 
Teff(E) 



(29) 



This formula gives us the escape time from the nonequi- 
librium potential barrier with the effective temperature 
which depends on the energy of the particle. 

In Fig. m we show the escape time ([29]) computed for 
various values of the asymmetry coefficient g as a func- 
tion of the applied voltage V. For all values of g the es- 
cape time decreases with increasing applied voltage. We 
also see that moving the molecular orbital closer to the 
Fermi energy of the electrodes facilitates the dissociation. 
When the noise is balanced by the viscosity, i.e. when 
the fluctuation-dissipation relations is forced for the nu- 
clear degrees of freedom, -?^ the dependence of the reac- 
tion rates on asymmetry coefficient is more complicated. 
For example, it was shown that the chemical reaction 
can be catalysed or stopped depending on the direction 
of the electric current. "'^^ In contrast, the absence of the 
fluctuation-dissipation relation leads to a rise of the effec- 
tive temperature Toff (2;) (or Tos^E)) which always over- 
rides the increase of the nonequilibrium potential barrier. 
We also see from the figure that for asymmetric coupling 
the escape time can be controlled by changing the applied 
voltage polarity, i.e., by replacing g hy 1/g. 



IV. CONCLUSIONS 

In this paper, we have formulated and solved Kramers 
problem for current-induced, out of equilibrium chemical 
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FIG. 6: Escape time computed for various values of the asym- 
metry coefficient g as a function of the applied voltage. The 
left and right panels correspond to different positions of the 
effective molecular orbital energy relative to the electrode self- 
energy eoff(a;min) ~ —0.36 + Vg. All other parameters are the 
same as those in Fig. (2] 



reactions. As a specific example, we have considered the 
dissociation of a H2 molecules induced by the tunnelling 
current from a STM tip. We have proposed a model 
Hamiltonian for the system and, by projecting out the 
electronic degrees of freedom, we have derived a Langevin 
equation for the time evolution of the reaction coordi- 
nate. The Langevin equation leads to Fokker-Planck dy- 
namics with an effective temperature which depends on 
the reaction coordinate. The reaction rate for the dissoci- 
ation is computed by solving the Fokker-Planck equation 
for the escape time. In the resonant tunnelling regime, 
two processes play equally important roles in the molec- 
ular dissociation: the decrease of the potential energy 
barrier and reaction-coordinate-dependent local heating 
of the molecule. The current induced catalytic reduction 
of the barrier accompanied by the local heating is an al- 
ternative mechanism for current-induced bond breaking, 
which is very different from the widely accepted paradigm 
of merely pumping vibrational degrees of freedom. 



Appendix A: Parametrization of the Hamiltonian 

In our calculations we take e{x) in the form of the 
bonding orbital for the H2 molecule in the gas phase^^ 



e{x) 



1 J(x) + K{x) 

2 ^ 



1 + Six) 

where the functions J{x), K{x), S{x) are given by 



(Al) 



J{x) — e 
S{x) = e 



'2x 



1 + - 1 K{x) = -e-'^il + x), 

X I X 



1+X + 



(A2) 



The Coulomb interaction between electrons Uq {x) is cho- 
sen such to reproduce, within the equilibrium {n{x) = 
1, Vg =Q) mean- field Hamiltonian ([2]), the energy of the 
H2 molecule obtained within the molecular orbital theory 
(see Appendix on page 543 in Ref. fiol). i.e.. 



{H-^l) = Emo{x) = 2e{x) + Uc{x) 



(A3) 



Here we neglect the nuclear kinetic energy contribution. 
Therefore, Uc{x) can be written as follows 



Uc{x) = 



^ + 0.5Jp(x) + Kp{x) + 2L{x) 
(l + 5(x))2 



(A4) 



where the functions L{x), Jp{x), Kpix) are listed in Ta- 
ble 10.8 of Ref. m 



Appendix B: Derivation of the Langevin equation 
for the reaction coordinate via NEGF 



The Langevin equation for nuclear degrees of freedom 
coupled to out-of-equilibrium electrons has recently been 
given in various contextaiii^Sr— . To make this paper self- 
contained, we provide here an explicit derivation which is 
directly adapted to the Hamiltonian under consideration. 
Our starting point is the Heisenberg equation of motion 
for the reaction coordinate as obtained from the mean- 
field Hamiltonian ([5]), 



dV, 



(Bl) 
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Here, the right-hand side contains the operator describ- 
ing the current-induced forces. Given that the molecular 
dissociation dynamics is slow compared to the electronic 
degrees of freedom, we calculate these forces within a 
nonequilibrium adiabatic approximation. In this approx- 
imation, the reaction coordinate is taken as classical. The 
electronic state should be computed for a given trajectory 
X = x(t) while in turn, the electrons affect this trajectory 
through adiabatic reaction forces. 
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In the nonequilibrium adiabatic approximation, one 
averages the force appearing on the right-hand side of 
Eq. (jBl|) over times which are long on the scale of the 
electronic dynamics, but short on the scale of the disso- 
ciation dynamics. The average adiabatic reaction force 
is thus given by the expectation value {e'^ff{x)a^^aa-)x(t), 
evaluated for a given trajectory x(t) of the reaction coor- 
dinate. The fluctuation-dissipation theorem implies that 
we need to include a fluctuating Langevin force 6f{t) in 
addition to the friction force. Hence Eq. (|B1I) becomes 



dVr 



N 



dx 



(B2) 



where we have introduced the lesser Green's function 



(B3) 



and the factor of two accounts for spin. Thus, except 
for the stochastic noise force, the current-induced forces 
are encoded in 2ie'gQ{x)Q^ . Below we demonstrate that 
the current-induced forces can be represented as a sum 
of two contributions: the adiabatic force and a velocity- 
dependent contribution. 



2ieMg< ^F{x)-ax)x. 



(B4) 



The variance of the stochastic force 5f{t) is gov- 
erned by the symmetrized fluctuations of the operator 
e[,ff (x)aJ.aCT. Given that the electronic fluctuations hap- 
pen on short time scales, 5f{t) is locally correlated in 
time. 



{5f{t)5f{t'))^D{x)5{t-t'). 



(B5) 



Since we are dealing with a mean-field Hamiltonian, D{x) 
can be evaluated using Wick's theorem, 



{5f{t)5f{t')) = 2[eMYg> it, t')g< {t\ t) , (B6) 



where 



g>{t,t') = -i{a,{t)al{t'))x^t) 



(B7) 



is the greater Green's function. These expressions show 
that we need to evaluate the electronic Green's function 
for a given classical trajectory x{t). 

We start with the Keldysh equation for lesser Green's 
function: 

= j d^i j dt2 g''{tM)^''{tiM)Q''{t2,t'), 

(B8) 

where the retarded and advanced Green's functions are 
given by the standard expressions 

g'^it, t') = ~id{t - t'){{a,{t), at , (B9) 



g''{t.t')^g'\t\t) 



(BIO) 



The Keldysh equation (jBSp involves the lesser self-energy 

i]<H = 2z ^/,(w)r,H, (Bii) 

a 

where fa{i^) = [1 + e'^°'^"~^°)]~^ is the Fermi-Dirac elec- 
tron distribution in the left (a = L) or the right (a = R) 
electrode, and Ta{ijj) is the imaginary part of the retarded 
self-energy due to interaction with the electrodes 



A„(cj) - ir„(w). (B12) 



The adiabatic expansion of Keldysh equation is carried 
out in the Wigner representation, given by 



'A{t + Tl2,t-T/2) 



(B13) 



for a general function A{t,t'), in which fast and slow 
time scales are easily identifiable. For the Green's func- 
tion g^'^'^ , the slow mechanical motion implies that 
^2) varies slowly with the central time t — 
, but oscillates fast with the relative time t = 
ti — t2- As usual, the Wigner transform of a convolu- 
tion C{ti,t2) = / dis ^(^1,^3)5(^3,^2) is given by 



C 



exp 
AB 



AB 



-d.AdtB - -dtAd.B, 
2 2 



(B14) 



where we have dropped higher order derivatives in the 
last line, exploiting the slow variation with the central 
time t. 

Expanding Eq. (jBSp up to the leading adiabatic cor- 
rection according to Eq. (jB14l) and taking into account 
that depends only on uj and is independent of the 
central time, we obtain g^ to first order in i, 

g< =G< +'-xe'^f,{x) [G<d^G> +G>d^G<] (B15) 

Here g denotes full Green's functions in Wigner rep- 
resentation, while G denotes the strictly adiabatic (or 
"frozen") Green's functions that are evaluated for a 
fixed value of a;: G^{x,uj) = [uj — ecff{x) — , 
S^M = EoS^(u;), G^ = (G^)t, G< = G«E<G^, 
G> = G< + G^ - G^. 

Let us now compute the current-induced forces appear- 
ing in the Langevin equation In the strictly adiabatic 
limit, i.e., retaining only the first term on the RHS of 
Eq. (|B15I) . g"^ ~ G^, we obtain the mean force 



F{x) = ~2s',s{x)n{x) 



(B16) 



where n{x) is the nonequilibrium population of the elec- 
tronic level 



n{x)^ [ ^G<{x,uj). (B17) 
J 2Tn 



The leading order correction in Eq. (jB15l) gives a 
velocity-dependent contribution to the current induced 
forces, which determines the friction ^(a;) in the Langevin 
equation. After integration by parts we obtain the ex- 
plicit expression 

^{x) = 2[eMY j ^G<{x,u^)d^G>{x,u). (B18) 

We evaluate the noise intensity D{x) (jB6|) for the stochas- 
tic force 5f to the lowest order in the adiabatic expansion, 
so that 

D{x)^2[eMY j '^G>{x,u^)G<{x,u^). (B19) 

Expressions Eqs, (IB18IB19I) can be simplified if we use 
the wide-band approximation for the electrodes. Namely, 
the wide-band limit employs that the retarded self-energy 
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is energy independent, (w) = —iTa- In that case 

G<[x,u:)^2^- (B20) 

and we we obtain the following expressions for the friction 

C(x) = 4r / E (-5-/a)ra (B21) 

a=L,R 

and the noise intensity 

/d(jj 
^Q(w,x) r„r„'/a(i-/a')- 
aa'=L,R 

(B22) 

Here, the factor Q{lo,x) is given by Eq. (|13l) . 
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